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Abstract. 

We present a cosmological model in which the metric allows for an inhomogeneous 
Universe with no intrinsic symmetries (Stephani models), providing the ideal features 
to describe a fractal distribution of matter. Constraints on the metric functions are 
derived using the expansion and redshift relations and allowing for scaling number 
counts, as expected in a fractal set. The main characteristics of such a cosmological 
model are discussed. 

PACS numbers: 98.80.-k, 98.65.Dx, 0.5.45.Df 
1. Introduction 

One of the most interesting results of observational cosmology over the last two decades 
has been the discovery of large scale structures and voids in matter distribution. 
To date, however, a clear definition of an homogeneity scale in matter distribution is still 
lacking. Moreover, this controversy concerns the morphological features of the galaxy 
clustering, which has been found by several authors to have fractal characteristics at 
least up to a certain scale (Sylos Labini, Montuori & Pietronero 1998; Wu et al. 1999). In 
this case, a suitable cosmological model is lacking, given that the Cosmological Principle 
cannot be applied. 

Many authors have attempted to develop alternative models, following two different 
approaches. The first relies on a perturbative scheme, through the superposition of a 
fractal density perturbation to the FRW background, both in Newtonian gravity (Joyce 
et al. 2000) and in general relativity (e.g. Mittal & Lohiya 2001). The second concerns 
a more general relativistic study of the metric, used to match a fractal number counts 
relation, and extensive analytical (e.g. Humphreys et al. 1998) and numerical (Ribeiro 
1992) computations have been performed. 

In any case, for both approaches, only Lemaitre-Tolman-Bondi (LTB) inhomogeneous 
metric has been considered, i.e. spherically symmetric solutions to Einstein field 
equations. 
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Even if this is certainly a very interesting study, it is not the most suitable tool to study 
a fractal distribution. 

As a matter of fact, in a fractal set of points the number density decreases on average 
from any point of the set. In a spherically symmetric distribution around a point, this 
is true only if you measure the density centering on that point. 

In this paper, the above restriction is relaxed and an inhomogeneous metric with no 
intrinsic symmetries is presented. 

The starting metric is recovered from the class of Stephani solutions (Stephani 1967). 
Specifically, the aim of this work is to extend the analytical characterization of a fractal 
matter distribution to a Stephani Universe, whose geometrical (e.g Krasinski 1983) and 
thermodynamical (e.g. Quevedo & Sussmann 1995) description has already been widely 
discussed and will not be pursued here. 

In Section 2 the general properties of the Stephani model are recalled while the 
consequences on the redshift expression and expansion are presented in Section 3. In 
Section 4 we derive the fractal constraints on the metric functions and the main results 
are discussed in Section 5. 

Throughout the paper we have assumed natural units (G = c = 1) and the metric 
signature (+,-,-,-). The Greek indices refer to the whole components, while the Latin 
letters are just for the spatial part. 

2. The model 

The metric of the Stephani Universe can be written as follows (Krasinski 1997): 

ds ? = D 2 df 2 _ &^l( dx 2 + dy 2 + d2 ?) (1) 

where : 

V = 1 + h(t){[x - x„(t)f + [y- m (t)f + [z- z (t)\ 2 } (2) 

D = F(t) (*-*)= F^ V - (3) 
KJ \V R) VdtR y ! 

k(t) = \C\t) - l/F\t)]R\t) (4) 

C, F, R, xq, yo, z are arbitrary functions of time and dots denote time derivative . 
Eq.(4) is equivalent to the Friedmann equation, although it must be stressed that the 
analogy is only formal and dictated by the same derivation (through Einstein field 
equations), but the function R(t) is not the equivalent of the scale factor a(t) in FRW 
models. 

As can be seen from eq.(l) and more precisely from the expression for V in eq.(2), the 
model is inhomogeneous and has no intrinsic symmetries, unless x ,y ,z are strictly 
constant. In this case xo,yo,zo become the center of spherical symmetry, even if the 
metric is still not homogeneous. 

Moreover, the distinction between open and closed Universe is more subtle than in FRW; 
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indeed k(t) is a function of time and so it is not fixed as in FRW model. However, it is 
worth remembering that fixing k(t) is not a general requirement of Einstein's theory of 
gravitation. 

Let us consider a perfect fluid; the stress-energy tensor is: 

T^u = (p + p)u^u v - pg^. (5) 

The Stephani metric (Eq.(l)) fulfills Einstein field equations if the matter fluid has no 
shear (a = ^a^a^ = 0), no rotation (uo = ^uj^uj^ = 0) and mass is flowing along 
t-coordinate lines with: 

u» = ^ (6) 

F RV — VR K ' 

«* = -£pD ti , u° = (7) 

* = = (8) 

where 5% is the Kronecker symbol. w M is the 4 — velocity, the notation H means the 
spatial derivative and ; the covariant derivative. In the FRW model, 9 = 3H where H 
is the Hubble constant. 

It is important to note that in the model the matter does not flow along geodesies, as 
follows from eq.(7); this is a significant difference with FRW model (which have u 1 = 0) 
(Krasinski 1983). 

Turning to thermodynamics properties, Einstein field equations with a perfect fluid 
source, the conservation law {Tf^ = 0) and the contracted Bianchi identities coupled to 
the matter conservation equation (nu^-^ = 0) can be reconciled in an equation of state, 
if the thermodynamical quantities are defined as follows: 

8vrp(t) = 3C 2 {t) (9) 

8Mt) = -3C 2 (t) + 2C\t)C^ - I) (10) 
3(p + p) 

n 4/3 W 

S = S + a,a=-^ (12) 
satisfying the Gibbs-Duhem equation (Quevedo & Sussman 1995): 



dS = ^ 
T 



d ^— ^ + p d 



n 



(13) 



As underlined by other authors (Krasinski, Quevedo & Sussman 1997), the physical 
meaning of the equation of state restricts the freedom in choosing the metric functions, 
though a defined correspondence among their properties and the physical quantities 
cannot be set, in general. However, since this work is mainly focused on a fractal model 
characterization, we will assume that the thermodynamic scheme holds and the metric 
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constraints will be derived independently. 

Again, a deep difference with respect to FRW models emerges, namely the equation of 
state is not barotropic, as it is evident from Eq.(9) and Eq.(lO). 

The fact that the equation of state depends on position, as indicated by eq.(9)-(12), 
and the absence of any intrinsic symmetry make the model particularly suited in the 
description of a fractal distribution of matter. 



3. The redshift equation and the expansion flow 

By definition, the redshift Z is determined along a photon path (written in terms of the 
affine parameter A) as follows (Ellis et al. 1985): 

(u»k u ) ( dt\ ,(dt\ 

(i+ ^ s ^=UL / uL (14) 

On the other hand, the definition of redshift refers to the ratio of the time shift along 
observer's world line with respect to the shift along a general world line, namely: 

+ (15) 

A more formal manner to find the aforementioned relation is given by considering the 
photon wave vector k^ 1 , which satisfies: 

dr^ 

fc% = (17) 

therefore, with no loss of generality, normalizing with respect to the observer's value, it 
yields: 

(1 + Z > = = L D = JW^V-R (18 » 
It is worth noting that the above redshift equation is dependent from the position, which 
is another important difference with respect to FRW models. 

In addition, other authors (Ribeiro 1992) have noted that more information about 

photon motion can be gained by evaluating the Lagrangian of the metric in eq.(l) 

and then solving for the Euler-Lagrange equations: 

d dL dL . dx^ , , 

, x» = — . (19) 



d\ dx^ dx^ d\ 
Inserting the past light cone expression relating the derivatives of the time and spatial 
parts of the metric (as follows by setting ds 2 = 0) into the Euler-Lagrange equation for 
the time component, a few algebra leads to the time motion of photons: 



ft_ 

d\ 2 



+ 



DDR V 



D D 3 RD 2 VD 2 




2 



TV =0. (20) 
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In the motion of a photon the affine parameter A is equivalent to the proper time r, so 
that, from eq.(15) and eq.(18), eq.(20) gives the redshift constraint on the metric: 

or equivalently: 

Comparing the V,R,F dependence in the redshift equation (eq.18) with the expansion 
expression in eq.(8) and including the redshift constraint (eq.(21)), yields: 

l = -jrh)^ 1+z) - (23) 

The left-hand side of the above equation is the equivalent of the Hubble parameter in 
the FRW models (3H — 9). In this case, the corresponding relation is: 

E ™ s -a-jrh)i (l + z) (24) 

where a(t) is the scale factor. 

Though a close correspondence between the two expression is apparent, they are quite 

different. In the case of an inhomogeneous model, the expansion is linked to the matter 

distribution since (1 + Z) (eq.(18)) depends on the position via the function V. 

In other words, different regions of the Universe can undergo different expansion 

histories. 

This is evident from the expansion law, which appears to contain an additional term. 
Indeed, the general redshift-distance relation is (Ehlers 1993): 

Z = + a^e^ 51 - u^5 ± x^ (25) 

where e M is the unit vector joining the observer and the emitting source, whose distance 
is 51, and directed towards the source. The position vector 5±x^ is defined in the 
following way. Let be the tensor projecting the tangent vector-space at each point 
perpendicularly onto the three dimensional subspace orthogonal to u\ namely: 

K = 5» + u»u v . (26) 

Considering a displacement 5x u , we can define the position vector as: 

5 ± x» = h»5x u (27) 

and then: 

51 = (g, v 5 ± x»5 ± xn 1/2 - (28) 

The shear contribution vanishes in the present model, therefore, by using the redshift 
constraint and introducing the acceleration in terms of the metric functions (eq.(7)), 
eq.(25) can be written as: 

Z=^- 2 5l + ^5x\ (29) 
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where ^-5x l is the kinematic acceleration evaluated at the position fixed by 51. 
The first term is the analogous to FRW Hubble term (albeit it depends on the position), 
since -j^ = H and SI is the comoving line element (see next section, eq.(40)). On the 
other hand, the additional second term is a new dipolar one, that is purely cosmological, 
i.e. not due to peculiar velocities. 

A few consequences of this will be further analyzed in the last section, nevertheless, it is 
worth noting that the present model is suitable to describe a fractal distribution, since 
all the cosmological parameters depend on the position. 

Using eq.(18) and eq.(29) , the metric function D must satisfy the following differential 
equation: 

^ + ^-1 + 1 = (30) 

and its solution can provide the space-time expression for the redshift. 
Eq.(30) basically describes a nonlinear wave motion with a damping factor. For instance, 
if the one dimensional equivalent is considered, eq.(30) has the general shape of the 
so-called Burger's equation, which is an approximate description of one dimensional 
turbulence, with a damping term. 



4. Fractal constraints 



The density of sources along the 51 direction after an affine parameter displacement (dX) 
at some point P is (Ribeiro 1992): 

d 2 N = 5l 2 dn [n(-k»u fl )] p d\ (31) 

where 51$ is the observer area distance (51 2 , = jj^) and n is the density of radiating 
sources in a subtended angle dflo. 

The density of sources can be found through the proper luminous matter density, 
assuming that the basic structures are galaxies with nearly the same mass Mq\ 

n = £=- = . (32) 

Using eq.(31) and eq.(32), a comoving observer measures a number of sources 
corresponding to: 

where integration is performed over the time At spent by the photon along its path; 
on the other hand, as a more physical interpretation, it directly links the motion of 
the photon to the intervening matter distribution, as the expression in eq.(9) for C 2 (t) 
shows. 

This value can be matched to the expected fractal value in the following way. 
Averaging over a sphere centered at the observer location, a volume and a volume density 
can be defined: 

V(di) = ^irdf (34) 
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PM) = -^j- (35) 

where this time the luminosity distance di corresponding to 5l was introduced, since it 
is the observational distance measure. It can be related to the observer area distance 
by means of the redshift factor (1 + Z) through the following: 

di = ' (1 + Zf = 5l (l + z f- (36) 

A fractal distribution of matter is characterized by a power-law scalingfj]: 

N(d l )=N [-±) (37) 




pM = ^if[r h ) ^ 

where lh is a transition scale corresponding to homogeneity and df = 3 — 7 is the fractal 
dimension. 

From the above equations and using the measured number counts, it yields: 

1-7 

(1 + Z ) 2{3 -^ (39) 



which links the unknown metric functions with observable quantities Nq, Mq, 7, lh and 
is the fractal constraint on the metric. According to eq.(9), the metric function C(t) is 
related to density evolution p(t), so that the physical meaning of the fractal constraint 
is to relate redshift and matter distribution in the proper way over each t = const 
hypersurface at each position. 

Evaluation of the metric terms at the homogeneity scale (7 7 (/h, Z(// l ))) simplifies the 
above relations and involves the knowledge of the fractal parameters (7, lh), the number 
of galaxies N and the typical galactic mass (M^), which are provided by the analysis 
of the galaxy catalogues, and the redshift factor, whose expression can be inserted from 
the solution of eq.(30). 

Otherwise, the value of Sl must be obtained using the proper-distance aperture of the 
subtended area, which follows from the line element for the metric in eq.(l). The line 
element along the i-th axis at a fixed timed to is therefore: 

^ w -rvS^ <4o) 

and can be found as (Krasinski 1983): 

2i?(t )[£;(to)n 1/2 arctan Q[A;(to)] 1/2 ^ ) k(t Q ) > (41) 
R{t )Sxi k{t ) = (42) 

f This scaling law should strictly apply only in an Euclidean space-time, since no general relativistic 
extension to it has been derived, so far. However, it has gained observational support, so it is used here 
without any theoretical speculations. 
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In addition, the value of N can be obtained by setting the matching condition to a FRW 
metric as lh is approached (Humphreys et al. 1998), i.e. by equating the expression for 
N(5l) evaluated at lh with the FRW analogous expression at the corresponding radial 
comoving homogeneity scale r^ (Coles & Lucchin 1995): 

, r , , , , f r h n[t(r)}a[t(r)]r 2 

N FRW (r h ,k FRW )=4n \ WJ 1 '' dr (44) 

Jo (1 - k FRW r 2 ) 

in terms of the FRW scale factor a(t) and curvature constant k FR w = — 1,0, +1. 
The above relation can be expanded and restated using the number counts continuity 
(na 3 = n ao), if the number of particles in the lapse time t — t(r) is kept fixed (i.e. no 
particle creation/destruction and evolution), as follows (Coles & Lucchin 1995): 

( r h 1 \ 
N FR w(r h , k FRW ) ^ Aim al f — - — k FRW 5l 5 + O J . (45) 

The condition: 

N FRW (r h , k FRW ) = N(l h , k(t)) (46) 



provides: 



N = Annual - ^k FRW r h 5 + OJ . (47) 



It should be emphasized that eq.(46) and eq.(47) allow a comparison between the basic 
FRW parameters (do, k FR w) and the ones involved in the present model, as they should 
represent the large scale structure in an equivalent way at the homogeneity scale. 
Although the fractal constraints provides a link among metric functions, it is more 
appealing to infer an analogous expression depending on the luminosity distance, being 
it the proper observational distance measure. 

In this case, the same route can be followed, but an important modification at the very 
beginning must be set. As discussed by other authors (Celerier & Thieberger 2001), if 
the observer area S(d a ) expressed in terms of the angular distance d a and the area S(di) 
resulting from the luminosity distance di are matched, then, using eq.(36), it yields: 

S(d l )=S(d a )(l + Z)\ (48) 

Again, considering eq.(36) and the solid angles ft(d a ), Q(di) respectively corresponding 
to S(d a ), S(di), by differentiating the above equation, an aberration effect on the solid 
angle is found: 



n{di) = n(d a ) 



-4/ -r 

JAt 1 



6/3 



-dt 



Q(d a )A( 9 -,z\ (49) 



I At (1 + Z) 

and the aberration function A (|, Z^j introduced, as can be seen, describes the effect of 
the expansion on the subtended angle, i.e. the integration of expansion in proper time 
provides the aberration of the solid angle. 

Using such a condition, the number counts in eq.(33) can be expressed in terms of the 
luminosity distance. In general, the latter is not a well defined quantity when referring 
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to cross sectional area perpendicular to light rays (SIq in eq.(31). The actual value of 
N(di) can be obtained by inserting eq.(36) and eq.(49) in eq.(31), then: 

JV(di) = H^r I ° 2{t \ dt. (50) 

By defining: 

we get the new fractal constraint: 

which indeed is now dependent on luminosity distance, redshift and expansion, namely 
the observational quantities. 



5. Discussion and conclusions 



The motivation for this work has been to describe a cosmological inhomogeneous 
model, which could represent a fractal distribution. Since spherical symmetry does 
not completely describe a fractal distribution, we have considered Stephani models, i.e. 
inhomogeneous models with no intrinsic symmetries. 
Our results are the following. 

First, the redshift expression depends on the position; of course, this character is present 
in any related cosmological parameter, such as the expansion and the Q value. 
In the FRW, expansion is given by the Hubble parameter, which is the only term 
governing the Hubble law. On the converse, in the present model, the Hubble law has an 
additional dipole term due to acceleration, which vanishes in the standard homogeneous 
cosmology (the quadrupole term due to shear vanishes as in FRW). Usually, the observed 
dipole term is ascribed to peculiar velocities; in the present model would add a pure 
cosmological contribution to it. 

The existence of a kinematic acceleration in inhomogeneous models has been invoked 
by other authors (Pascual-Sanchez 2000) as the origin of accelerated expansion of the 
Universe, as found by recent SNae la data (Perlmutter et al. 1999; Riess et al. 1999). 
This could be considered as a viable alternative to the proposed interpretation within 
FRW scheme, namely a positive cosmological constant or vacuum energy or quintessence. 
Therefore, the striking implication of the coupling of space and time dependences in 
redshift is that different regions of the Universe could undergo different expansion 
histories. 

Actually, it should be stressed that the present model would lead to a deeply new 
interpretation of observations and many cosmological and astrophysical contexts should 
be reconsidered. 

A second result is that we have explicitly introduced a fractal characterization of the 
model. 
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The density depends on the redshift, which depends on the position. In addition, it 
fulfills the scaling law imposed by the fractal distribution. Then, the density changes 
according to the fractal spatial scaling law and the redshift. This approach allows to 
map the density evolution on t = const hypersurfaces, applying the fractal constraint. 
The homogeneity scale lh plays an important role in the model. We recall that it is 
the scale at which the Stephani model should match the FRW one, i.e. the scale at 
which homogeneity and isotropy of the Universe are recovered. The greatest information 
supply should come by estimating the parameters at that threshold length and imposing 
the FRW/fractal matching condition. 

Summarizing the whole set of constraints, the present model can explicitly express D, 
F and C from functions in the metric. 

The key information still missing concerns probably the most extreme assumption of 
the model k(t), i.e. curvature time evolution, that is no more kept fixed and bound 
to -1, 0, +1 values as in FRW Universe. Once k(t) is specified, R and V can be fully 
determined and a self-consistent cosmology can be built. 

So far, we have listed the main differences between the FRW model and the one 
proposed here. However, an important point against inhomogeneous models is their 
incompatibility with high degree of isotropy of the Cosmic Microwave Background 
(CMB). Indeed, it was demonstrated that an isotropic thermal radiation in an expanding 
Universe implies a FRW cosmology (EGS theorem, Ehlers, Green & Sachs 1968). The 
point is that the EGS theorem was demonstrated for geodetic and non-accelerating 
observers, which is not the case in Stephani models. If acceleration is not vanishing, 
an inhomogeneous cosmology allows a thermal radiation with high degree of isotropy 
as observed (Ferrando, Morales & Portilla 1992) and specific Stephani solutions can be 
built without any violation to the CMB restriction (Clarkson & Barret 1999). 
Nevertheless, it is to stress that limitations to the degree of inhomogeneity of the model 
must be compatible with, for instance, CMB isotropy and Hubble law local trend; for 
such a purpose, the derived redshift and fractal constraints must be tuned. 
Actually, a more subtle and rather philosophical argument comes into play when facing 
with any inhomogeneous cosmology, that is the reliability and interpretation of the 
Cosmological Principle. 

Moreover, the present model poses serious questions about the Copernican Principle, as 
well, which states that all observers are equivalent in the Universe. This is the basis, 
for instance, of the EGS theorem and is far from being straightforwardly fit by a fractal 
structure. 

A last comment concerns the observations, most of which can be reproduced by FRW 
cosmology. As a conservative statement, it could be noted that each Stephani solution 
includes FRW as a limiting case, so that observations cannot in general rule them out. 
The motivation for keeping FRW as a description of the Universe is actually that they are 
the easiest theoretical model. This led to prefer any modification to come from outside 
the cosmological framework (e.g. inflation, dark matter) and not in the cosmology itself. 
Therefore, the theoretical choice has been overwhelming any observational basis, whose 
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actual interpretation is not a priori excluded and could be reanalyzed within the class 
of inhomogeneous models. 
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